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ABSTRACT 

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James 
Glerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 

160, 250, 350, 450, 500 and 850 jam. The disk is resolved at 70, 160 and 450 jam. In addition to the 
planetesimal belt, we detect thermal emission from AU Mic’s halo for the first time. In contrast to 
the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the 
disk is 3.9 x 10“^ and its mm-grain dust mass is 0.01 (±20%). We create a simple spatial model 

that reconciles the disk SED as a blackbody of 53 ± 2 K (a composite of 39 and 50 K components) and 
the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is 
consistent with the “birth ring” model explored in earlier works, i.e., an edge-on dust belt extending 
from 8.8-40 AU, but with an additional halo component with an surface density profile extending 
to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force 
to blow out grains. Eor stellar mass loss rates of lO-lOOx solar, compact (zero porosity) grains can 
only be removed if they are very small; consistently with previous work, if the porosity is 0.9, then 
grains approaching 0.1 jam can be removed via corpuscular forces (i.e., the stellar wind). 

Subject headings: (stars:) circumstellar matter, individual (AU Mic) 
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1. INTRODUCTION 

Debris disks are one of the most prevalent signposts 
that a stellar system succeeded in building up planetary 
scale bodies during the protoplanetary disk phase. De¬ 
bris disks are collisionally sustained distributions of plan- 
etesimals, smaller rocky bodies, and dust around main 
sequence (and more evolved) stars. Because dust grains 
can be removed from the system through various physi¬ 
cal processes, their presence is direct evidence of an un¬ 
seen population of larger planetesimals, and potentially 
planets, in orbit around the star. At current sensitivity 
levels, debris disks are found around > 20% of nearby 
solar and A-type stars dEiroa et al.||2Q13| |Thureau et al. 
20141 |Matthews et al.||2014b| IVIattnews et al. in prep), 
with enhance d detection rates arou nd younger stars, e.g., 
the Pleiades: G orlova et al. (|20 06 ) fin d 25% for B and A- 
type stars while' bierchio et al. '(^'OlO) find 32 % for solar 
types: the 1 0-20 JVIyr old bco^U en association: |Chen et ^ 


(|2011|) and Chen et al.| (|2012|) find rates from 25 — 33%, 
varying with spectral type ; and ensembles of A star pop¬ 


ulations: Su et al. (2006) find a rate of 32% at 24 and 
70 jam. There is, howev er, a relative paucity of disks de¬ 
tected around M stars (Matthews et al.||2014b). Recent 
surveys for debris disks have revealed very low detection 
rates for M star hosted debris disks compared to earlier 
type stars ( [Low et ^|2005 Gautier et al.||20 07). For ex- 

th~ 


ample, data from the DEBRlb survey with the Herschel 
Space Observator'^^ieYesl just t wo disks in a sample of 
89 observed M stars within 8.6 pc (iMatthews et al.|2015, 
in prep). The mass sensitivity to M star disks in existing 
surveys, however, has not y et matched those of ey lier 
spectral types (see Fig. 2 of Matthews et al. 2014b). 

Several factors affect the detectability of disks around 
M stars. For example, grain removal by stellar winds 
is mo re efficient around M stars (e.g., jPlavchan et al. 
2005), suggesting that small grains (i.e., < 1 jam) may 
be removed from these systems at a r ate higher than ex¬ 
pected purely from radiation forces ( [Augerean fc Beust 


2006 Matthews et al. 2007). This ettect may explain 
why the highest detection rate of 13lg% is found for 
a combination of submillimeter studies since these are 


sensitive to l arger grains (Lestrade et al. 2006). 


For- 


brich et al. (2008) find that tor the 30 — 40 Myr old 


cluster NGC 2547, the detection rates of M star disks 
at 24 jam exceeds that of G and K stars of the same 
age, suggesting that at least around very young stars, M 
star disks may be just as detectable as disks around ear¬ 
lier type stars, cons istent with the bright disk detected 
around the M3IV e dPecaut fc Mamaj^ 2013) star AU 
Mic, the M3IVe (jPecaut ^ M amajek||2013|) star TWA 
7’s disk ([ Matthews et al.||2007|) , and the recently identi¬ 
fied disk ([K ennedy et al.^Oldt ar ound the M4 star Fo- 
malhaut C (|Mamajek et al. 5013). Nevertheless, a disk 
of high fractional luminosity has also been detected and 
resolved arou nd the significantly o lder, multiple planet 
host GJ 581 ([Lestrade et al.[j2012[). Therefore, AU Mic, 
as the first and youngest nearby M star to have a detected 
disk, retains particular importance as a representative of 
its class. 

AU Mic (GJ 803, HD 197481) was the only M star 

Herschel is an ESA space observatory with science instru¬ 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 


detected to have a far-IR excess above its stellar pho¬ 
tosphere (indicative of c ircumstellar dust) wi th IRAS 
(Faint Source Catalogue, Moshir & et al.[[1990[). A sub- 


millimeter excess was detected with photometric-mode 
observations with the SCUBA camera, and the emis¬ 
sion was no t resolved in a relatively shallow mapping 
observation (Liu et al. 2QQ 4[). An excess wa s also seen 
with the CSO at 350 jam ^Ghen et al. 2005). Ground- 


based imaging by Kalas et al. (|2004p and subsequent 
high-resolution im aging and po l arime try f rom the Hubble 


Space Telescope by Krist et al. (2005) and Graham et al. 
(2007) revealed it to be the second debris disk spatially 


resolved at optical wavelengths, hosting an edg e-on disk 
that extends to over 100 AU in radius. The IGrahaml 


et al. (2007) HST polarization study revealed evidence 
of a change in the polarization properties at a radius of 
^ 35 A U and a dearth of mi cron-sized grains interior to 
40 AU. Wilner et al. (2012) imaged the disk with the 
SubMillimeter Array (I8MA) and most recently, imaging 
with the Atacama Large Millimeter/submillimeter Array 
(ALMA) yielded a well-resolved disk image modeled as 
a narrow “ birth ring” or “parent b elt” of planetesimals 


at 40 AU ( [MacGregor et ah 


2013[). The ALMA data 
do suggest a reasonably wide dust belt extending inward 
to ^ 9 AU though the dust surface density is strongly 
peaked near 40 AU and the inner edg e is poorly con¬ 
straine d. In recent STIS observations, [Schneider et al.[ 
(2014) report an “out-of-plane bump” on one side of the 
disk at ^ 13 AU, interpreting this as a dust density en¬ 
hancement, in contrast to the ALMA data that revealed 


no significant asymmetries. Schneider et al. (2014) also 
observe a brightness asymmetry between the two sides of 
the scattered light disk interior to 15-20 AU, a region that 
had not been cleanly imaged in earlier scattered light de¬ 
tections. Finally, the ALMA data also revealed an unre¬ 
solved excess at the position of the star which MacGregor 
et al. note could be attributable to unresolved emission 
from an aste roid-like warm belt near the star. ICranmerl 


et al. (2013), however, suggest coronal thermal heating 
could alternatively account for the observed excess at the 
position of the star. 

AU Mic attains additional importance as a member of 
the Beta Pictoris Moving Group since, as such, it has 


([Binks & Jeffries 

2014 Malo et al. 

2014 Mamajek & 

Bell 

2014), revised upward from 12^4 Myr (Zuckerman 


et a 


pc. 


van Leeuwen 


2007), AU Mi c is a favored target for 


study. As described in detail by Wilner et al. (2012), its 


disk also shares many qualitative characteristics with the 
P Pictoris disk, including its ed jge-on geometry an d ex¬ 
tended scattered light emission (Kalas et al.[[2004[). The 
age of AU Mic is the currentl y favored epoch for the for¬ 
mation of terrestrial planets ( Chambers | p014 Raymond 


et al. 2012). While M stars do host planetar y systems, in¬ 


cludin g fne mu ltiple planets arou nd GJ 876 (Ri vera et al. 
2005[), GJ 581 ([Ud ry'e t al.|2007 ) and GJ 676A ( [Anglada- 


Kscude Tuom 2012| ). No planets have been detected 


around AU Mic, although Schneider et al. (2014) do sug¬ 
gest a potential planetary origin for the observed asym¬ 
metry in scattered light. 

AU Mic is similar to other bright debris disks in that 
it possesses a bright “halo”, seen in scattered light im- 












































































































































































































AU Mic: Herschel and JCMT resolved imaging 
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ages. In debris disk systems where the smallest grains 
are blown out of the system as soon as they are created, 
and larger grains are unaffected, there necessarily exists 
an intermediate size range where newly created particles 
are placed on eccentric orbits. The specific sizes of these 
particles depends on the stellar luminosity and mass loss 
rate, but in general they are smaller than 10 /im. These 
particles have pericenters within the parent belt (or birth 
ring), but apocenters extending to the maximum allowed 
by the local environment. Thus, they form a small-grain 
“halo” that surrounds the parent belt. In the case of AU 
Mic, the radiation force is relatively weak and the halo 
is created by the radial force exerted on the dust by the 
stellar wind, which is thought to be 10-100 times greater 
than the Solar wind (|Strubbe & Chiang||2006t |Augereau 
& Beust||2006|). A key prediction is therefore that ha¬ 
los should be relatively faint at mm wavelengths where 
the small grains emit inefficiently, and thus only the par¬ 
ent belt is detected. This hypothesis has been confirmed 
by mm-wave observatio ns of AU Mic ( Wilner et al.|20lT 


MacGregor et al.|2013 ), as well as for other systems (e^ 


Vega, 13 Fie, and HR 8799; [Holland et al.|1998| |Williams 


& Andrews] 2006 Dent et al. 2014 respectively), thus 
explaining how debris disks can have different appar¬ 
ent sizes at different wavelengths. For both Vega and 
HR 8799 the halo has only been detected in mid/far-IR 
emission (iSib thorpe et al.|2010||Sn et al.|2009| [Matthews 
et al.||20l4a| r 


We present the first resolved far-infrared images of 
AU Mic’s debris disk from the Herschel Space Obser¬ 
vatory's PACS (Photodetector Array Camera and Spec- 
Poglitsch et al.|2010 ) instrument at 70 /am and 


trometer, 


160 jam and lower resolution submillimeter images from 
the SPIRE (Spectr al and Photometric Imaging REceiver, 
Griffin et al.||2010 ) camera at 250 /am, 350 /am and 500 


jam. I'hese data were taken during Guaranteed Time, as 
part of the Disk Evolution Key Program (PI: Olofsson). 
In addition, we present 850 /am and 450 /am data from 
the SCUBA-2 Observations of Nearby Stars (SONS) Sur- 
vey and a PI program , taken with the SCUBA-2 camera 
Holland et al.|[2Ql^ . We describe our observations and 
' ata reduction in § and present the data in § HI In 
§ we discuss the new constraints these data place on 
models of the disk emission, namely the temperature of 
the planetesimal belt and halo and the extent of the halo 
in thermal emission. We summarize our conclusions in § 

m 


2. OBSERVATIONS 

A summary of deta ils of the observations from Herschel 
(jPilbratt et al.|2010|) and the James Clerk Maxwell Tele- 
scope (JCIVIT') is shown in Table [h The relative sensitiv¬ 
ities of Spz^2;er/MIPS24, Hersd^/VKCS and SCUBA-2 
observations to the AU Mic disk in terms of fractional 
luminosity (Ldisk/A*) against blackbody radial dust lo¬ 
cation are shown in Eigure[2 


Blackbody temperature (K) 

490 160 49. 16. 4.9 



Figure 1. Relative sensitivities of SpitzerMlVS (24 /xm), Herschel 
PACS and SCUBA-2 to dust emission at the level of AU Mic’s disk, 
shown as a black dot. Any disk above each instrument’s sensitivity 
curve will be detectable to that instrument. The best fit temper¬ 
ature and fractional luminosity of AU Mic render undetectable to 
MIPS at 24 fim (as reported by|Chen et al.|2005|), but detectable 
to PACS at 70 and 160 jim as well as bCUBA-z at 450 and 850 

/j,m. 

total observing time was 3.8 hours. The images at each 
observed wavelength are shown in Eigurej^ 

The PACS observations used scan legs of 7.4' with a 
scan-leg spacing of 38". Medium scan rate maps contain¬ 
ing 11 scan legs were repeated 11 times in each scanning 
direction to achieve the required map depth. Data were 
obtained at both 70 /am and 160 /am simultaneously. 

The nominal cross-linked large map observing parame¬ 
ters were used for the SPIRE observations, with a smaller 
map of 4 arcmin by 4 arcmin being executed. In total, 
ten repeat maps were used to reach the required depth. 

The data were processed using version 13.0 of the Her - 
schel interactive pipeline environment (HIPE,|Ott|2010). 
The standard pipeline processing steps were used tor 
both the PACS and SPIRE data. Versions 69 and 13.1 of 
the PACS and SPIRE calibration products were applied 
respectively. The final maps have pixel scales of 1, 2, 6, 
10 and 14" at 70, 160, 250, 350 and 500 /am respectively. 

As part of the data processing, the PACS time-lines 
were first masked at the source location, as well as in 
other areas of bright emission, and then high-pass fil¬ 
tered to reduce the impact of 1/f noise. The filtered 
data were then converted to maps using the ’photPro- 
jecU task. Maps were likewise made from the SPIRE 
data using the ’naiveMapper’ task. No filtering was per¬ 
formed on the SPIRE data. 

Observations of a Tan and Neptune, adjusted to the 
correct spacecraft position angle, were used as the model 
PACS and SPIRE PSEs respectively, and used for model 
image convolution. 

2.2. JCMT Observations 


2.1. Herschel Observations 

Large area photometric mapping observations were 
performed using the PACS and SPIRE cameras at 70 
/im, 160 /im, 250 /im, 350 /am and 500 /am. The large 
scan-map mode was used and nearly perpendicular cross- 
linked scans were performed with both instruments. The 


Observational data presented in this paper were also 
taken using the SCUBA-2 camera (Holland et al. 2013) 
on the JCMT. The data were o btained both aspa rt of 
the SONS JCMT Legacy survey (Phillips et al. 2010) and 
a PI program. The observati ons were taken with t he con¬ 
stant speed DAISY pattern (|Holland et al.||2013|), which 
maximizes exposure time and provides unitorm cover- 
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Table 1 

Observations Log 


Obs. ID 

Observing date 

Mode 

Duration [s] 

PWV [mm] 

Herschel Space Observatory 

1342196038 

9 May 2010 

PacsPhoto 70/160 (scan 135) 

5478 

— 

1342196103 

9 May 2010 

PacsPhoto 70/160 (scan 45) 

5478 

- 

1342193786 

5 April 2010 

SpirePhoto 

2906 

- 

James Glerk Maxwell Telescope 

SGUBA-2 

22 April 2012 

daisy scan 850/450 

3910 

1.06 


23 April 2012 

daisy scan 850/450 

3947 

0.83 - 0.92 


17 May 2012 

daisy scan 850/450 

1960 

0.92 


7 June 2012 

daisy scan 850/450 

1895 

1.4 


16 August 2012 

daisy scan 850/450 

3787 

0.86 


19 August 2012 

daisy scan 850/450 

3790 

0.83 



30 20 10 0 -10 -20 30 20 10 0 -10 -20 30 20 10 0 -10 -20 30 20 10 0 -10 -20 -30 
Offset / arcsec Offset / arcsec Offset / arcsec Offset / arcsec 



Figure 2. Far-infrared and submillimeter maps of AU Mic from Herschel and SCUBA-2. North is up and east is to the left. Different 
surface brightness scales are used for each map, and the pixel scales used are 1" at 70 /rm, 2" at 160 /rm, and 6", 10" and 14" at 250 /rm, 
350 fim. and 500 /rm, respectively. The pixel scale of the JCMT maps is 1", and the maps are smoothed to a half- (850 fim) or full-width 
(450 fim) Gaussian. The green contours show the 3 — cr level in each of the maps. The 1 — cr rms levels are 0.9 mJy beam“^, 9.0 mJy 
beam“^, 5.8 mJy beam“^, 6.3 mJy beam“^, 6.0 mJy beam“^, 6.8 mJy beam“^ and 0.9 mJy beam“^ from 70 through 850 fim. The rms 
levels for the SPIRE data are the confusion limits of the instrument {Herschel Observers’ Manual). The background source, “BG”, is well 
isolated from the AU Mic disk emission at 70 jim and surrounded by a 3 — cr contour at 70 /j,m and 160 jim (labeled). 


age in the central 3' diameter region of a field. The 
total integration time was just over 5 hours, split into 
10 separate ^ 30 minute observations. Observing con¬ 
ditions were generally excellent with precipitable water 
vapour levels less than 1 mm, corresponding to zenith 
sky opacities of around 1.0 and 0.2 at 450 and 850 mi¬ 
crons respectively (equivalent to JCMT weather “grade 
1”; i" 225 GHz of < 0.05). The data were calibrated in flux 
density against the primary calibrator Uranus and also 
secondary calibrators CRL 618 and CRL 26 88 from the 
JCMT calibrator list (Dempsey et al. 2013), with esti¬ 
mated calibration uncertainties amounting to 10% at 450 
jam and 5% at 850 jam. 

The SCUBA-2 data were reduced using the Dynamic 


Iterative Map-Maker within the STARLINK SMURF 
package (Chapin et a l.|[2M3) called from the ORAC-DR 
automated pipeline (|(Javanagh et al. 2008). The map 
maker used a configuration file optimized for known posi¬ 
tion, compact sources. It adopts the technique of “zero- 
masking” in which the map is constrained to a mean 
value of zero (in this case outside a radius of 60" from 
the center of the fi eld), for all but the final interation 
of the map maker (Chapin et al. 2013). The technique 
not only helps convergence in the iterative part of the 
map-making process but suppresses the large-scale rip¬ 
ples that can produce ringing artefacts. The data are 
also high-pass filtered at 1 Hz, corresponding to a spa¬ 
tial cut-off of ^ 150" for a typical DAISY scanning speed 
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RA offset (a rcsec) 


Figure 3. SCUBA-2 imaging of the AU Mic disk. North is up and 
east is to the left. Composite image showing the extension of the 
disk along the known disk orientation (shown by the orientation of 
the scale bar) at 450 p,m in greyscale with an overlay of the 850 
lim contours. The contour levels are RMS levels of 3-cr (white), 
6-a (green), 9-cr (light blue), 12-cr (blue) and 15-cr (black). 


of 155"/s. The filtering removes residual low-frequency 
(large spatial scale) noise and, along with the “zero- 
masking” technique, produces flat, unifor m final images 


large ly devoid of gradients and artefacts (Chapin et al 


2013). 

To account for the attenuation of the signal as a result 
of the time series filtering, the pipeline re-makes each 
map with a fake 10 Jy Gaussian added to the raw data, 
but offset from the nominal map center by 30" to avoid 
contamination with any detected source. The amplitude 
of the Gaussian in the output map gives the signal at¬ 
tenuation, and this correction is applied along with the 
flux conversion factor derived from the calibrator ob¬ 
servations. The final images were made by coadding 
the 10 maps using inverse-variance weighting, re-gridded 
with 1-arcsec pixels at both wavelengths. The final im¬ 
ages at both wavelengths have been smoothed with a 7" 
FWHM Gaussian to improve the signal-to-noise ratio. 
The FWHMs of the primary beam are 7.9" and 13.0" at 
450 jam and 850 /im, respectively. 


3. RESULTS 

Figure shows the observed thermal emission on the 
sky toward AU Mic at seven wavelengths. The maps 
have been cropped to show only a 1' x 1' field centered 
on the star. AU Mic is detected at all wavelengths. The 
disk is resolved at 70 and 160 jam with Herschel. Two- 
dimensional Gaussian fits to the JCMT data yield sizes 
of 1670 X 876 at 450 jam (at PA ^ 135°, aligned with 
the scattered light images) and 1679 x 1474 at 850 /im, 
which suggest that the source is resolved along the major 
axis at 450 jam (as also indicated by the difference in the 
peak and integrated flux densities in Table ^ but only 
marginally at 850 /am (see composite Figure^. 

The PACS flux densities were measured using aper¬ 
ture photometry, with radii of 17" at 70 /am and 34" at 
160 jam. We used aperture corrections of 0.81 and 0.85 
at 70 and 160 /im, which were derived from a large set 
of calibration observations processed in the same way as 
the data. At 160 /im, this aperture includes the back¬ 
ground source, whose point source flux was estimated in 



Figure 4. Optical image of AU Mic from HST showing the ori¬ 
entation and extent of the scattered light disk. North is up and 
east is to the left. The 70 /j,m emission from Herschel is overlaid. 
The disk is clearly resolved at 70 /xm, and the disk orientations are 
consistent, to the limits of the Herschel resolution. The yellow box 
shows the position of the 70 /j,m source detected with Herschel at 
the epoch of the HST observations. The position is close to, but 
not coincident with, the optically visible galaxy approximtely 2" 
to the southeast. 


the modelling described below and has been subtracted 
to estimate the AU Mic flux density alone. The RMS 
noise levels were estimated by integrating the flux within 
beam-sized circular apertures (^ 10 at 70 /im), spaced 
by the FWHM beam, over the central 2' diameter area. 
However, the PAGS flux densities have an uncertainty 
that is largely set by the instrumental calibration, with 
some additional uncertainty at 160 jam to allow for sub¬ 
traction of the background source. At 70 jam and 160 /im, 
the absolute flux calibration accuracies are 3% and 5% 


respe ctively {Herschel Observers’ Manual; Balog et al. 


2014). Due to the depth of the image, the measurement 


errors are negligible compared to the calibration errors, 
so the combined residual sum of squares uncertainties are 
the same as the calibration errors. We attempted to mea¬ 
sure SPIRE flux densities using PSF fits. Our flux ex¬ 
traction at SPIRE wavelengths, however, is severely im¬ 
pacted by our inability to separate the flux density of the 
AU Mic disk from nearby background objects, so we do 
not quote any fluxes. The JCMT flux density measure¬ 
ments were made using 20" radius apertures, with the 
same method for uncertainty estimation. At these wave¬ 
lengths, however, the calibration uncertainty (^ 10%) is 
relatively unimportant. 


Table reports the peak flux densities and integrated 
flux densities of AU Mic at each observing wavelength. 
As well, we list the PSE-fit flux density of the “BG” 


background source identified in Eigure These mea¬ 
surements are found to be in good agreement with those 
derived from the image model of AU Mic described be¬ 
low. 


Rebull et al. (2008) report flux densities of 143 ±2 mJy 
at 2 4 /im, 205 ± 8 mJy at 70 /im, and 168 ± 20 at 160 


jam. Plavchan et al. (2009) report updated flux densities 
of 155.2 ± 3.2 mJy and 223 ± 26 mJy at 24 and 70 jam. 
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Figure 5. SED of AU Mic and its disk. The fit to the stellar 
photosphere is given in blue; the disk fit is shown in red, and the 
composite stellar + disk spectrum is shown in black. Black symbols 
are observed flux densities and grey symbols show star-subtracted 
values. Inverted triangles show upper limits (the 100 p,m value is 
from IRAS). In the absence of resolved imaging, a single component 
disk would be a reasonable interpretation of this SED, for which 
the excess emission is well characterized by a pure blackbody. 

The photospheric flux at 24 jam is 150 ±2 mJy so there is 
no significant excess at this wavelength. Our measured 
flux density at 70 /im from Herschel is consistent with the 
measured Spitzer flux density, while we measure a higher 
flux density (but consistent within 2-cr) than Spitzer at 
160 jam. 

The only significant source of emission within 1' of AU 
Mic is located at PA = 244.0° and separation ^ 25". To 
determine if this is a background galaxy, we checked the 
HST F814W observations made with ACS/WFC in 2004 
(GO-10228; PI Kalas). Figure 4 shows the 70 jam Her¬ 
schel contours overlaid on the HST image. No galaxy is 
detected in these optical data that corresponds to the 70 
jam source to the southwest of AU Mic (position at epoch 
of optical data is given by the yellow box), although the 
galaxy in the optical image is quite close (^ 2", the point¬ 
ing accuracy of Herschel). The absence of a coincident 
galaxy is not unexpected given the relative sensitivity to 
infrared faint galaxies of Herschel compared to optical 
observations. 

The flux density distribution (which we loosely call the 
spectral energy distribution, or SED) of AU Mic and its 
disk is shown in Figure To derive the stellar spectrum 
we compiled UBV, Hipparcos/Tycho-2, 2MASS, Spitzer, 
AKARI, and WISE photometry up to 12 jam. We ex¬ 
cluded the U — B color to avoid potential issues with vari¬ 
ability due to flaring at the shortest wavelengths. The 


best fit to the B T-Settl stellar atmosphere models (Al- 
lard et al.||2012|) yields Te// = 3600 ± 20 K, a radius of 
0.83 and a luminosity of 0.1 Lq. Low-res IRS Spitzer 
data range from 5 — 14 jam and are consiste nt wi th the 




stellar photosphere, but not shown on Eigure 
We then subtracted the photospheric model from the 
flux densities at longer wavelengths, and fitted a pure 
blackbody to this star-subtracted photometry. The re¬ 
sulting model is shown in Eigure where the disk has a 
temperature of 53 ± 2 K and a Actional luminosity of 

A high-resolution IRS Spitzer scan was examined but was not 
useable. 


3.5 X 10“^. That the disk spectrum is well fit by a pure 
blackbody is surprising for two reasons; first, most debris 
disks have an emission spectrum s teeper than Rayleigh- 
Jeans at lo ng wavelengths (e.g. 


et al. 


Lm|[2QQ4 


2012), and second, the disk is known to be 


2008 Caspar 


extended from scattered lig ht imaging dKalas et al. 12004 


IKrist et al.||2005l IMetchev et al. 


itzger- 


ald et al. 2007t|Graham et al.|2007||Schneider et al.|26l4p 

so must comprise dust at a wide range of temperatures. 

The goal of the spatial modelling in the next section is 
to use a simple model to reconcile the disk extent and 
blackbody spectrum. 

3.1. Modeling 

We use a simple model to interpret the Herschel^ 
SCUBA-2, and ALMA observations. Preliminary tests 
find that the Herschel 70 jam image is significantly more 
extended than the parent belt imaged with ALMA at 1.3 
mm, so the basic requirement of the model is that it rec¬ 
onciles the different extent of these images and produces 
a blackbody-like spectrum. As is already known, the so¬ 
lution is that the greater radial extent at 70 jam originates 
in the halo of small grains also seen in scattered light, and 
that these are not seen with ALMA because these grains 
emit inefficiently at wavelengths significantly larger than 
their physical size. Thus, our model comprises two com¬ 
ponents; the first is the parent bodies for which we us e 
the ALMA modelling results ofjMacGregor et al.|(|2013|), 
and the second is a halo whose basic properties are to be 
informed by previous theory work and determined from 
the modelling. 

With only six photometric points in the spectrum, and 
a only few beams of resolution in the new Herschel and 
JCMT images, our modelling approach is physically mo¬ 
tivated but simple. It has been used previ ously to model 
many Herschel-ie^o\ved debris d isks (e.g., Kennedy et al. 
2012a|b Matthews et al. 2014a). Eor each disk compo- 
nent a single azimuthally symmetric 3D dust distribu¬ 
tion is used, which is simply a small scale height disk 
with a power-law surface density dependence between 
the inner and outer radii. At each radial location, the 
disk emission is assumed to arise from a modified black¬ 
body, with a power-law radial temperature dependence. 
This approach is therefore largely empirical, and the de¬ 
rived blackbody parameters can be subsequently com¬ 
pared with more detailed dust models to draw conclu¬ 
sions about the dust properties, in particular their typi¬ 
cal size. The limited resolution also precludes derivation 
of comparative radial profiles, as was done for HR 8799’s 


halo (Matthews et al. 


2014a). 


Practically, we generate high resolution images at each 
observed wavelength, as viewed from a specific direction 
to produce the desired disk geometry. We can include an 
arbitrary number of disk components, as well as point 
sources to model unresolved sources such as warm in¬ 
ner disk regions or background sources. Each image is 
then convolved with the appropriate point spread func¬ 
tion (PSE), which may be an observation of a bright cali¬ 
bration source (i.e., PAGS, SPIRE), or simply a Gaussian 
with specified width(s) and position angle (i.e., JGMT, 
ALMA), and then resampled to the resolution of the ob¬ 
servation. We model the ALMA image and not the vis¬ 
ibilities, which is acceptable given the good uv coverage 
and dynamic range of the ALMA data. The models are 
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Table 2 

Measured fluxes 


Wavelength 

[fim] 

Peak Flux 
[mjy beam“^] 

Integrated Flux 
[mJy] 

Fbg 

[mJy] 

P'phot 

[mJy] 

Disk Flux 
[mJy] 

70 

91.6 ±2.7 

231 ±7 

11.6 ±0.35 

19.6 

219 ±7 

160 

176 ±8.8 

228 ± 15 

47.1 ±2.3 

3.6 

226 ± 15 

450 

32.3 ±5.5 

49.2 ±8.5 

non-detection 

0.44 

35.4 ±8.5 

850 

12.9 ±0.9 

14.0 ± 1.5 

non-detection 

0.10 

12.5 ±1.5 


Note. — The columns contain Q) the observing wavelength; (2) the peak flux intensity as 
measured from the maps of Figure!^ the uncertainty represents the 1 —cr RMS in the maps; (3) 
the integrated intensity measured^ aperture photometry; (4) the flux of the SW background 
source as derived from a PSF fit, color-corrected; (5) the predicted flux density of the stellar 
photosphere, color-corrected; and (6) the disk flux, color-corrected. The uncertainties on 
the extrapolated photospheric fit are of the order of a few percent and therefore negligible 
compared with the measured uncertainties. 


fit to the data using a combination of by-eye variation of 
parameters and least-squares minimization. Since we do 
not search all parameter space, our best-fit models are 
not necessarily unique, but must be considered a reason¬ 
able interpretation. 

Our model is guided by previous modelling work. The 
parent belt is modeled using the structure derived by 
MacGregor et al. (2013) based on ALMA data, here us- 
ing a surface density distribution between 8.8 AU and 40 
AU with a radial power law dependence of H oc r^-^. The 
inner edge is not tightly constrained. MacGregor et al. 
find that the inner edge of the parent belt could be as far 
out as 21 AU, sugg esting a much narrowe r disk width, 
though we note that Schneider et al. (2014) find evidence 
of a stello-symmetric warp in the disk with an outer edge 
of 15 AU, which supports the idea that some component 
of the disk has an inner edge closer than 21 AU. Note 
that we are calling this component the parent belt for 
convenience, and that the extent inward of 40 AU may 
be due to stellar wind drag from a true source region that 
is much narrower, but likely still concentrated at ^ 40 
AU. The emission properties of grains in this belt are as¬ 
sumed to be blackbody-like, and follow Tpb = To,pb/\/f, 
meaning that the su rface brightness profile m atches the 


power-law found by MacGregor et al. (2013). Because 


radiation and wind forces are relatively weak for AU 
Mic, small dust can reside in the planetesimal belt and 
To,pb can be higher than expected for blackbodies, so we 
leave it as a free parameter. We also add the unresolved 
point source component at the stellar position seen in the 
ALMA data ([MacGregor et al.|2013|), though this source 
is too faint to be detected at any other wavelength (see 


3.2). 


To reproduce the Herschel 70 /im image we include 
the halo componen t, with a surface density p rofile fixed 


to S oc (c.g., Strubbe & Ghiang 2006|), which ex¬ 

tends from 40 AU to 140 AU. ('Lhe outer edge is poorly 
constrained by the data and fixed at 140 AU.) The sur¬ 
face density of this component is forced to join smoothly 
to the planetesimal belt, and the grain properties are 
allowed to vary via their temperatures and emission 
spectra. Their spectra are modified blackbodies (i.e., a 
Planck function multiplied by (Aq/A)^^ for A > Aq), but 
we reduce the number of parameters by enforcing a tem¬ 
perature law suited for small grains, Thai© = 7o,haio^~^^^- 
We vary Aq, /3a, and To^haio as free parameters. As noted 
above we expect Aq ^ 100 /im and /3 a > 0, so that the 
halo is detected at 70 /im, but not at 1.3 mm with ALMA. 


These values are consisten t with measured values of An 
and 0 for ot her disks (i.e.. Booth et al.||2013 Matthews 
et al.||2014a ). 


finally, we include the point source visible to the west 
of AU Mic in Figure This source is included in part 
to ensure it does not bias the model at 160 /im, but is 
also used to derive a flux density to subtract from the 
160 jam aperture flux derived above. We also include 
the central point source seen in the 1.3 mm ALMA im¬ 
age. With this simple model we have four main param¬ 
eters; one for the parent belt temperature and three for 
the halo component temperature and spectrum. There 
are in addition eight RA/Dec offsets for the Herschel^ 
JGMT, and ALMA images, but these are relatively well 
constrained. The aim of the model fitting is therefore 
essentially to find the relative weight of the two compo¬ 
nents at each wavelength, thereby empirically deriving a 
coarse spectrum for each in a way that is independent of 
assumptions about grain properties. We found that the 
relatively poor resolution (excepting the ALMA data) 
means that the parameters are poorly constrained due 
to degeneracies, in particular for the halo component. 
However, all solutions we found are sufficiently similar 
that the ultimate interpretation is the same. 

Figure ^ shows the original images, the parent belt and 
halo model components, and the model-subtracted resid¬ 
uals (for just the planetesimal belt and the full model) 
at wavelengths of 70 /im, 160 /im, 450 jam and 1300 jam. 
First, to make it clear that the halo is detected, the fourth 
column shows the data with only the parent belt sub¬ 
tracted. In this column the brightness of the parent belt 
is 1.8 and 1.4 times brighter than in the final model at 
70 and 160 jam. The 70 jam residuals clearly show that 
the disk is more extended than the parent belt, and that 
a component with greater extent (i.e., the halo) is re¬ 
quired. The fifth (rightmost) column shows the residuals 
with the addition of this component, and that our parent 
belt -1- halo model reproduces the data well. Thus, the 
halo contributes significant emission at 70 /im, but little 
at longer wavelengths. 

At all wavelengths, few residuals above the 3-<j level 
are seen. The negative residuals at 160 jam appear to lie 
along the scan directions (at roughly ±45° relative to N) 
suggesting that they are artefacts of the imaging. The 
positive residual to the SE lies beyond the disk extent and 
is not seen on the opposite side of the star, so should not 
be considered as a deficiency of the model. The ALMA 
residual map has a bright point source offset from the 




























Figure 6. Resolved images of the AU Mic disk at (top to bottom) 70 p,m, 160 p,m, 450 p,m and 1300 p,m (left column, as in Figure 
compared to the planetesimal belt (second column) and halo (third column) components. The fourth column shows residuals (and ±2- ana 
3-cr contours) when only the planetary belt model (multiplied by 1.8, 1.4, 1 and 1 for 70, 160, 450 and 1300 /rm, respectively) is removed 
from the data and the firth column shows the residual when the full (parent belt + halo) model is subtracted from the data, and includes 
additional background components, which are most visible in the 160 p,m image. The color scale for the first three panels along each row 
is the same to show the relative contribution of each component. The color scale on each row is different, scaled to near the peak ffux in 
the observed image at that wavelength. The color scale of the residual images are scaled to the same Ter level. While the residuals at 450 
p,m are suggestively symmetric, they are not aligned with the known disk orientation and are barely at the 2-a level. In addition, a point 
source at the stellar position was included in the model at 1.3 mm, consistent with the unresolved excess reported by [MacGregor et al.| 
( |2013| >. The residual to the south of the star at 1.3 mm is assumed to be an unrelated background source. 


ste Har position. This is a b ackground source identified 
by MacGregor et al. (2013|) as well. It is located 0''3 
west and 178 south of the star and is detected with a 
flux density of 0.15 mJy (a 5-cr detection). 

The main derived model parameters are To,pb = 245 
K, To,halo = Ibl K, Ao = 12 /im, and /3x = 1. Because the 
structure of the parent belt is fixed, To,pb is moderately 
well constrained by the contribution required at 70 jam. 
The halo properties, which are our primary interest, are 
less well constrained because is it only strongly detected 
at 70 fim. However, the conclusion that Aq is shorter 
than ^70 /im, and that > 1 is robust because the 
halo is required by the images to have little contribution 


at 160 yum and beyond. Given the conclusions of previ¬ 
ous studies that the halo is populated by small dust on 
eccentric orbits these parameters are consistent with our 
expectations. The specific value of Aq = 12 /im suggests 
^1 /im sized grains dominate the emission, but a poor 
constraint on this parameter means that the grains could 
be larger or smaller, though by no more than an order of 
magnitude. We return to the likely size of bl owou t grains 
when considering small-grain dynamics in 


4.2 


Figure shows the star-subtracted SED and the disk 
model, and here the relative contributions of the two 
components as a function of wavelength can again be 
seen. We note that the ALMA fluxes (1300 /im) lie some- 
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the disk, we would not be able to detect a coronal ther¬ 
mal heating contribution to the total flux in the Herschel 
data. 


3.3. Mass of the Disk 

A direct estimate of the mass of dust in AU Mic’s debris 
disk can be made from its submillimeter flux densities. 
Debris disks are optically thin at these wavelengths, and 
so the mass of the disk is directly proportional to the 
emission, following the relation: 


Md = 


F, d? 

Ku Bi,{Td) 


( 1 ) 


In the Rayleigh-Jeans limit, equationreduces to 


= 5.8 X 10-1° 


Fi,[mJy\ (d[pc])^ (A[/xm])^ 
K^[cm?g-^\ Td[K] 


( 2 ) 


Figure 7. Flux distribution of the disk around AU Mic. Subtrac¬ 
tion of the stellar fit from Figure reveals the distribution of the 
disk itself. Mid-infrared upper limits are shown from WISE (22 
/jm), Spitzer (24 pun) and IRAS (25 iim). Far-infrared detections 
at 60 yum (IRAS) and 70 /jm and 160 fim (Spitzer and Herschel) 
are shown. Submillimeter detections are from the CSO (350 p-m), 
JCMT (450 and 850 p,m) and ALMA (1.3 mm). The model com¬ 
ponents of the planetesimal belt (“parent belt”) and the halo are 
shown. 

what above the fit, but are consistent with the model at 
the 2 — (7 le vel. We discuss the spectral index of the 
disk in § |4.3| below. As noted above, at 70 /rm the emis¬ 
sion from the planetesimal belt and halo components is 
similar, but at longer wavelengths the planetesimal belt 
dominates. The planetesimal belt is dominated by emis¬ 
sion at 40 AU, where the grains have a temperature of 39 
K. The temperature of the halo component at that sepa¬ 
ration from the star is slightly higher (50 K), but poorly 
constrained due to the lack of resolution and degeneracy 
with Ao and px. The fact that these two components 
yield a composite spectrum that is a simple blackbody is 
a reminder of the power of resolved imaging, since much 
information can be hidden within a single component 
SED. 


3.2. Central Asteroid Belt? 


MacGregor et al. (2013) found evidence of an excess 
above the photosphere at the position of the star in their 
high resolution ALMA imaging, and we include this com¬ 
ponent in the ALMA model image. The residuals from 
fitting our dust model to the Herschel and JCMT images 
in Figure 6 show no sign of any other unresolved excess 
above the stellar photosphere, so the best we can do is 
place an upper limit on the emission from such a cen¬ 
tralized component of the p oint-source sensitivity a t a 
given wavelength (i.e., 3 — a, MacGregor et al.||2013). It 
is also possible, given the considerably lower resolution 
of our observations, that such emission, if present, has 
been incorporated into the belt emission and removed, 
which adds unc ertainty to the flux de nsity of the belt. 

From Fig. 4 of Cranmer et al. (2013), the coronal model 
of AU Mic predicts a flux of a tew /ijy around 100 jam 
with a relatively flat distribution shortward of 1300 jam 
where the ALMA detection was made. Our point source 
sensitivity is on the order of 1 mJy at 70 fim and ~ 10 
mJy at 160 /im, meaning that, even in the absence of 


While the AU Mic SED supports a single tempera¬ 
ture fit of 53 K, our image modeling requires a separate 
planetesimal belt of 39 K and a warmer halo, which con¬ 
tributes negligible flux density to the total emission at 
850 /im (see Fig. 7). This temperatur e is comparable to 
the 40 K adopteoDy Liu et al. (2004|) in the absence of 
an SED temperature fit. T'he color corrected flux density 
measured with SGUBA -2 at 850 /am (12 .5 mJy) is within 
the 1 — cr limit of the value measured by Liu et al. (2004) 
with SGUBA (14.4 ± 1.8 mJy), although the SCUBA 
measurement relied on single bolometer “photometry” 
mode, rather than a mapping technique, which could 
have missed some emission. 

Taking 12.5 mJy as the 850 /am flux density of the 
disk yields a slight ly lower 1 -mm y ain dust mass of 0.01 
M 0 relative to the Liu et al. (2004) measurement, for the 
same dust opacity of 1.7 cm^ g“-^. The dominant sources 
of uncertainty in the mass are the flux density and the 
temperature of the planetesimal component of the disk, 
amounting to a 20 % uncertainty in the disk mass. 

4. DISCUSSION 

We now briefly consider the dust properties in more 
detail, using realistic grain models to compare the prob¬ 
able dust temperatures and sizes with those derived from 
the modelling. 

4.1. Temperature of AU Mic Dust Grains 

The SED fit to the disk around AU Mic requires only 
a single temperature component. The emission is well 
fit (x^ is minimized) from mid-IR through submillimeter 
wavelengths by a pure blackbody with a temperature of 
53 ± 2 K. Figure shows, however, that the halo and 
planetesimal belt components could still be segregated 
in temperature, with the parent bodies of the belt colde r 
(^ 39 K) than the halo (see also Fitzgerald et al. ]2007l) . 

Figure shows the distribution ot temperature with 
grain size for three grain compositions at the location 
of the planetesimal belt. The temperature dis tributions 


were calcula ted as in Augereau et al. (1999) and Wyatt & 
Dent ( 2002 ). We used three difterent compositions: i) a 


mix of 1 part amorphous astronomical silicates to 2 parts 
organics, ii) pure astro nomical silicates, and iii) pure ice 
(Li & Greenberg 1997). The models shown in Figure]^ 
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Figure 8. Temperature of grains as a function of diameter for 
three dust grain compositions. For the planetesimal belt tempera¬ 
ture of 39 K, pure icy grains are ruled out for all considered grain 
sizes for compact, spherical grains, although other compositional 
mixes including ices could potentially satisfy the modeled temper¬ 
ature. See the text for an explanation of the grain compositions. 

are not porous. Porous grains do not have the peak in 
temperature near 1 /im, instead decreasing steadily to 
20 — 25 K from maximum temperatures similar to the 
non-porous case. It is clear that pure water ice grains 
are ruled out at 39 K for all grain sizes, although we 
note that other compositional mixes could be modeled 
to satisfy the temperature distribution we have derived. 


4.2. Small-grain dynamics 

Orbiting dust is subject to forces from both radiation 
and stellar winds. Around early-type stars radiation 
dominates, while stellar winds are thought to be more 
important for late-type stars such as AU Mic. These 
forces are commonly split into radial and tangential com¬ 
ponents. The radial component is called ’’radiation pres¬ 
sure” or the ’’radiation force”, and effectively reduces the 
stellar mass seen by the particle, with an analogous effect 
for the stellar wind. The tangential component is called 
Poynting-Robertson drag, and causes dust grains to spi¬ 
ral into the star, again with an an alogous (and much 
stronger) effect for stellar wind. See Burns et al. (1979) 
for a detailed review. 

The small grain halo seen previously in scattered light 
and confirmed here in the far-IR is similar in morphol- 
ogy to those seen ar o und other debris disks (e.g., (3 Pic: 
gantm_et a^ JT9^ ); Angereau et al. (| 2QQ1 ); HR 8799: 


Sn et al. ( |2U09^1Vlatthews et al. 


terns, the halo is attributed to 


). in other sys- 
le efiect of the radial 


component of the radiation force from the star, which 
increases the eccentricity of smallest bound grains. The 
key measure of this effect is the parameter the ratio 
of the radiation force to the gravitational attraction of 
the star (distinct from the mm-wave spectral slope Px 
used above). Grains with /3 > 1 are unbound, and grains 
liberated from a parent body on a circular orbit are un¬ 
bound if /d > 0.5. Figureshows /d as a function of grain 
diameter, calculated for amorphous silicates + organics, 
and ice as described above. The solid curves show (3 for 
the effect of radiation force on solid grains (left panel) 
and porous grains (right panel). It is clear that a range 
of /3 is possible due to the poorly constrained grain com- 


positio ns. Also, the assumption of Strubbe fc Chian^ 
(2006) that the radiation force blowout grain size exists 
(i.e., that P > 0.5 for any size) is not well founded, par¬ 
ticularly if the grains are porous. Essentially, the low 
luminosity of an M type star does not produce a strong 
enough radiation force to significantly affect the small 
grains. Even taking into account the flaring of the star 
is not enough to remove grains via the radiation force 
( Augerean fc Benst]|2006 ). 

Unlike early-type s tars, an M star can hav e a relatively 


strong stellar wind (Plavchan et al. 
from this particle wind causes a simr 


2005). The force 
ar enect as the ra¬ 


diation force and is abl e to reproduce the scattere d light 
halo around AU Mic ^ugereau &: Bens^|2006 ). The 
dotted lines in Eigure demonstrate the eftective 
parame ter for the stellar wind fo rce (using the prescrip¬ 
tion of Strubbe & Chiang 2006). The left panel shows 
that for compact grains the wind blowout size is very 
small, <0 .1 jUm for the favo red AU Mic mass loss rate of 
10 X Mq ( Cully et al.|l994 based on EUV flaring). Such 


small grains are extremely poor scatterers of light, even 
at optical wavelengths, so they would appear very faint. 


Angereau & Beust (2006) get around this problem by in- 
cluding flares, which increase the time-averaged mass loss 
rate to 300 x Mq . Such a high mass loss rate would how¬ 
ever make the system transport-dominated rather than 
collision-dominated, i.e., stellar wind drag would also 
become important, filling the region interior to 30 AU 
with small dust grains in a manner contrary to th e scat 


tered light obse rvations (Strubbe & Chiang 2006). Lim 


& White (1996) also point out that radio flares would be 


obscured if the stellar winds of M stars were many orders 
of magnitude more massive than the Sun’s. 

The right panel of Eigureshows that porosity can also 
play an important role. By increasing the grain porosity, 
the effect of wind force increases, in turn increasing the 
blowout grain size. Eor p = 0.9 the blowout size is an or¬ 
der of magnitude larger for the same mass loss rate, and 
the scattering efficiency of the grains greater, thus pro¬ 
viding the likely solution. Eor a stellar wind rate of ^100 
Mq the size of halo grains is roughly 0.1 to 1 pm. These 
grains have emission properties consistent with those de¬ 
rived from the image modelling ab ove. This conclusion o f 
high porosity compares well with Graham et al. (2007), 


Eitzgerald et al. (2007) and bhen et al. (|2009p who all 
find that porosity is necessary to explain the scat tering 
and polariza tion properties of the grains, although |Shen| 


et al. (2009) find that a composition of random aggre¬ 
gates require a por osity of just ^ 0.6 co mpared to the 
0.9-0.94 required by|Grahani et al.|(|2007|) for Mie theory 
applie d to spherical grams or agy egates. 

The |Strubbe & Chiang| (|2006|) model of the disk as 
a collision-dominated, narrow birth ring from which the 
smallest grains are blown out does a good job of explain¬ 
ing the halo and lack of small grains interior to 30 AU . 
The recent ALMA results of MacGregor et al. (2013), 
however, show that there are indeed larger grains interior 
to the presumed birth ring. Given the 1300 pm observ¬ 
ing wavelength, those observations should be dominated 
by grains with a ^ Xl27r, or 200 pm. Drag forces, which 
can act on larger grains over long timescal es, are there¬ 
fore a likely cause of this interior emission. [Angereau fc 


Beust (2006) concluded that Poynting-Robertson drag is 
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Figure 9. Effects of Porosity. The ratio of radiative to gravitational forces (/3) as a function of grain composition is shown. For zero 
porosity grains, there is no grain size that reaches the blowout condition (/3 = 0.5) around AU Mic, although non-zero porosity can either 
decrease or enhance the /3 ratio, depending on grain size. For the higher porosity grains, the /3 values ultimately converge to be similar, 
following an inverse relation to grain size. Stellar wind effects on /3sw for I? 10? 10^ cind 10^ X Mq are shown as dotted lines. Empirical 
evidence for other stars suggests that AU Mic may have a stellar wind level as high as, but not much higher than, 10 X Mq (Strubbe & 
|Chiang|2006| >. 


negligible for AU Mic and that stel lar wind drag may 
evacu ate the inner regions (see also Strubbe & Chiang 
2006). However, more recent work suggests that stel¬ 
lar wind drag can fill the interior regions to some degree, 
without seriously violat ing the scattered light constraints 
dSchiippler et aL||2015|). Alternatively, the radially in- 
creasing surface density of this disk is reminiscent of a 
“self-stirre d” disk that collisionally grinds down from the 
inside out ([Kenyon fc Bromley[|2002 Kennedy & Wyatt 


2010). Schiippler et al.| (|^015p also hnd that the self- 
stirred model is plausible tor AU Mic’s disk. 

4.3. Spectral Index at millimeter wavelengths 

The data compiled in this work provide a long lever 
arm to measure the spectral index, a where oc 
is the emission from submillimeter through millimeter 
wavelengths. As noted above, the 1300 /im fluxes from 
the SMA and ALMA both lie above the nominal “best 
fit” line of the models produced by fits to the SED and 
the images. Comparisons of the 350/1300 spectral index 
and the 450/1300 spectral index yields a values ranging 
from ^ 1.5 to close to 2.0, which was the inferred spec tral 
index of 350 vs. 1300 (SMA) by Wilner et al. (2012). It 
is clear from the SED that these wavelengths are m the 
Rayleigh-Jeans part of the spectru m, which wasn ’ t defi¬ 


nite based on the data available to Wilner et al. (2012) 


and Caspar et al. (2012). 

A spectral index of ^ 2 is lower than most other de¬ 
bris d isks. Based on the compiled data of [Caspar et al. 
(|2012|), the three A stars in that study Pictoris, Vega 
and bomalhaut) have the highest value of a, which could 
suggest that the size distribution is different around A 
stars, being more co nsistent with a collisional cascade. 


Caspar et al. (2012) show that for a collisional quasi 
steady-state at a single temperature, the Rayleigh-Jeans 
tail of the SED should yield a spectral index of 2.65, for 
a self-similar collisional differential size distribu tion in¬ 
dex c) f 3.65. The shallower slopes measured by [Wilner 
et al. (2012) and a comparison of the 450 /im flux den- 
sity from SCUBA-2 of Table (which actually seems to 
lie somewhat below the SED nt) and the elevated flux 


densities measured by the SMA and ALMA at 1300 /im 
suggest a size distribution index of 3.0, i.e., that the size 
distribution is not the product of a collisional cascade. 
The shallow slope may however arise from the way that 
the parent belt and hal o components add to yield the 
total spectrum (Eig. 7, [Eitzgerald et al.[[2007|). Such a 
shallow distribution constraint has also been noted for 
two IR bright disks (HD 377 and HD 104680) not de¬ 
tected with the VLA at 9 mm (Creaves et al.[[2012). In 
AU Mic, it appears the grains rnay be very large (>> 1 
mm). Longer wavelength observations of the disk beyond 
1300 /im would be very beneficial to further constrain the 
size distribution. 

The spectral slope is also dependent on composition. 
Schiippler et al. (2015) present more detailed modeling of 
the AU IVlic disk, including investigation of the impact of 
compositional dependence on the derived spectral slope. 
They find a best fit for a combination of silicate, carbon 
and vacuum in equal measure, with little change in the 
spectral slope fit with small additions of ice or variations 
of these materials. Their test also provides clear evidence 
for porous grains. As in our analysis, their SED fits 
also underestimate the 1300 jam flux densities from SMA 
and ALMA, though as in our case, this result should 
be interpreted with caution since the exploration of the 
parameter space is limited. 

5. CONCLUSIONS 

AU Mic’s halo has been kno wn since the earlie st scat¬ 
tered light images of the disk ( Kalas et al.[[2004 ), which 
revealed a radial extent of 210 AU. ALIVIA imaging has 
co nfirmed a planetesim a l belt at ^ 40 AU, as predicted 
by St rubbe & Chiang (2006) and Augereau & Beust 
([2006[) and consistent with a break in the surface bright¬ 
ness profiles observed in scattered light (jMetchev et al. 


2005 Eitzgerald et al.[[2007 

Graham et aU 

2UU7(). Using 

Herschel and SCUBA-2 at t 

ne J CIVIL, we 

aave detected 


and modeled the thermal emission from both the plan¬ 
etesimal belt and the halo at wavelengths of 70 /im, 160 
/im, 450 /im an d 1300 /im, the latter be ing ALMA data 


as presented by MacGregor et al. (2013). 
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We present a simple spatial model that utilizes the ex¬ 
isting model of the planetesimal belt from ALMA imag¬ 
ing to reconcile the ostensibly single-temperature black- 
body SED of AU Mic’s disk with the presence of emission 
from the extended halo. The best fit model is consistent 
wit h the birth ring model ex plored in Wilner et ah (2012) 
and|MacGregor et ah (2013), a planetesimal belt extend- 
ing from 8.8-40 AU, but with the addition of a shallow 
surface density profile halo dominated by grains roughly 
1 fim in size. 

We observe no asymmetries in the disk images, and the 
residual images all show that there is negligible emission 
unaccounted for by a smooth disk. 

We confirm that AU Mic does not exert enough radi¬ 
ation force to blow out grains. We also find that for the 
inferred stellar mass loss rate of 10 times solar, compact 
(porosity = 0) grains can only be removed if they are very 
small; if the porosity reaches 0.9 or higher, then grains 
approaching 0.1 /im can be removed. This result sug¬ 
gests that a higher mass loss is favoured to place larger 
^1 /im grains in the halo, and a high degree of porosity 
in the grains of AU Mic, consistent with previous work 
on the scattering and polarization properties of the disk 
at optical wavelengths. 

The spectral index of the planetesimal belt of AU Mic 
may be more shallow than our modeling suggests, if the 
450 /im diminished flux density and 1300 /im elevated 
flux densities are real. The spectral index established 
from 350 /im and 850 /im data was already shallow at a 
value of ^ 2, but may be as low as 1.5, suggesting that the 
disk may have a grain size distribution inconsistent with 
that expected of a quasi steady-state collisional cascade. 
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